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Abstract. We establish locality estimates, known as Lieb-Robinson bounds, for the Toda lattice. In 
contrast to harmonic models, the Lieb-Robinson velocity for these systems do depend on the initial condition. 
■ Our results also apply to the entire Toda as well as the Kac-van Moerbeke hierarchy. Under suitable 

assumptions, our methods also yield a finite velocity for certain perturbations of these systems. 

o 

wo: 

^ . 1. Introduction 

Analyzing the dynamics of certain non-relativistic systems is crucial to understanding a number of im- 
£C) ' portant problems in statistical mechanics. For example, there is much interest in rigorously justifying the 

emergence of macroscopic non-equilibrium phenomena, like Fourier's Law or other forms of heat conduction, 
directly from a many-body Hamiltonian dynamics 0]. Future progress on fundamental questions such as 
these will require detailed information on the structure of the underlying dynamics. One of the goals of the 
Oh| present work is to investigate an approximate form of locality, often described by a Lieb-Robinson bound, 

^ fin for the dynamics corresponding to a class of integrable systems. 

Recognizing approximate forms of locality, despite the lack of a relativistic framework, has been essential 
in solving many intriguing open problems. One of the first mathematical formulations of a useful locality 
estimate was given by Lieb and Robinson in 1972 |24j . In this work, they demonstrated that the dynamics 
corresponding to quantum spin systems, with e.g. finite-range interactions, remains effectively confined to a 
£N) | "light" cone, up to corrections which decay at least exponentially away from the light cone. Recently there 

have been a number of improvements and generalizations of the original result [5TJ EQ1 EH QUI El ESI EH ES 
[391 EI] , and these new techniques have led to some interesting applications [171 H51 1521 1501 17| IH1 151 151 US ] . 
For a review of these results, we refer the interested reader to [331 ITS] . 
\ Shortly after the original result of Lieb and Robinson, it was shown in [26] , see also [9j [38] for more recent 

• ■ developments, that this notion of quasi-locality also applies to the dynamics of classical oscillator systems. 

Since it is this work that more closely pertains to the topic of the present article, we will discuss it briefly 
as follows. 

Consider a system of particles confined to a large but finite set Ac Zf To each site x G A associate 
a particle, or oscillator, with position q x € R and momentum p x 6 M. The state of the system in A is 
described by a sequence x = {(q x ,Px)}x&A, and the set of all such sequences, X\, is called phase space. A 
Hamiltonian, H, is a real- valued function on phase space. Given a Hamiltonian and a sequence x 6 X\ 1 
^ | Hamilton's equations of motion are: for each ieA, 

. . . 8H , . dH 

(1-1) QxW = ^— and p x [t) = - — 

opx oq x 

solved with initial condition {(q x (0),p x (0))} xe \ — x. For many Hamiltonians, this system of coupled 
differential equations can be solved for all time (and any initial condition). In this case, we denote by 
$ t , the Hamiltonian flow, i.e., the mapping that associates to initial conditions, the solution at time t: 
* t (x) = {(q x {t),p x (t))} xe A. 

Measurements of the system under consideration correspond to observables, where an observable A is a 
complex- valued function on phase space. For example, the position of the particle at site x € A corresponds 
to an observable Q x for which Q x (x) = q x . Let us denote by A\ the set of all observables over X\. Given 
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a Hamiltonian for which (jl.ip can be solved, the Hamiltonian dynamics at : A\ — > Aa is defined by 
at{A) = A o In this case, we can make time-dependent observations such as [at(Q x )](x) = q x (t). 
It is well known that the Hamiltonian dynamics is generated by the Poisson bracket, i.e., 

(1.2) j t a t (A) = a t ({A, H}) = {a t (A),H} , 

and here the Poisson bracket of two observables is the observable 

(13) { A B} = Y — ■ — - — ■ — 

teA ® qx ^ Px ® Px ^ qx 

In terms of this Poisson bracket, we can now describe the Lieb-Robinson bound. 

Let us fix finite sets X, Y C 1 d with X n Y = 0. Take A c 1 d finite, but large enough so that lUFcA. 
Consider two observables A,Be Aa with supports in X, Y respectively, i.e., e.g. A depends only on those 
q x and p x with It is clear that {A, B} — since A and B have disjoint supports. A Hamiltonian H 

satisfies a Lieb-Robinson bound if for some initial condition x there exist numbers fi, C, and v for which 

(1.4) \[{a t (A),B}](x)\ < c e -/*WW-«|t|) , 

where d(X, Y) is the distance between X and Y. In words, this estimate shows that for times t with 
v\t\ < d(X, Y) the support of a t (A) remains essentially disjoint from the support of B, up to exponentially 
small corrections in d(X, Y). The number v is often called the Lieb-Robinson velocity corresponding to H, 
and it represents a bound on the rate at which disturbances can propagate through the system. 

We will now briefly discuss what is known concerning Lieb-Robinson bounds for classical systems. The 
works in 26 , 9 , 38 prove Lieb-Robinson type bounds for a variety of harmonic systems, and it is shown 
that analogous results also hold for certain anharmonic perturbations. The strongest such result is in [5] 
and demonstrates that for an anharmonic model with a quartic on-site perturbation the relevant Poisson 
bracket decays to zero whenever the distance between the supports of the local observables grows faster than 
tlog a (t) for suitable a > 0. In 38], where explicit estimates on the Lieb-Robinson velocity were obtained, 
it is shown that the number v is independent of the initial condition x. Since the system is linear, this is not 
so surprising. 

Of course when it comes to anharmonic lattice oscillations one of the central objects is the famous Fermi- 
Pasta-Ulam-Tsingou (FPU) problem. It was first demonstrated by Zabusky and Kruskal [44] that one 
key ingredient to resolve the FPU paradox is the relation with solitons. Moreover, subsequently Toda 
presented an anharmonic lattice which possess soliton solutions and was later on shown to be integrable by 
Flaschka [HI Q3] and Manakov 125. . Clearly this naturally raises the question about Lieb-Robinson type 
locality bounds for the Toda lattice and the main goal of this work is to establish such bounds. Our estimates 
produce a Lieb-Robinson velocity that depends on the initial condition, see e.g. Theorem 12.31 We discuss 
this fact in the context of one-soliton solutions in Section [2T4l For a restricted class of initial conditions, we 
can prove the existence of a finite Lieb-Robinson velocity for a class of perturbations of the Toda lattice. 
We have two bounds of this type. The first, in Section |3~T1 establishes a result by directly mimicking the 
methods for the unperturbed system. The next, in Section l3~2l uses interpolation and shows that the velocity 
of the perturbed system can be estimated in terms of the velocity of the unperturbed system. Analogous 
results are shown to also hold for solutions of the Toda hierarchy, see Section [4] and Section [5] For a much 
larger class of initial conditions, we can prove a locality bound for perturbed systems, however, we do not 
establish the existence of a finite Lieb-Robinson velocity. This is discussed in Section [5] The final section, 
Appendix [Al describes a set of perturbations and corresponding initial conditions for which solutions remain 
globally bounded. This shows that the results obtained in sections 13.11 and [3.21 are not vacuous. 

2. Locality Estimates for the Toda Lattice 

In this section, we discuss the Toda Lattice and prove a Lieb-Robinson bound. We introduce the model 
in Section [2~T1 The crucial solution estimate, Theorem 12. 1[ is contained in Section |2~21 as well as some useful 
remarks. A Lieb-Robinson bound, see Theorem 12 .31 in terms of a large class of observables is proven in 
Section [2131 We end this section by comparing the velocity estimates in our Lieb-Robinson bound to known 
results for one-soliton solutions of the Toda Lattice. This is done in Section [2^41 
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2.1. The Toda Lattice. The Toda Lattice is a well-studied physical model and one of the prototypical 
discrete integrable wave equations. We refer to the monographs [11], [40], [43] or the review articles [22], [41] 
for further information. Like those discussed in the introduction, the Toda Lattice corresponds to a specific 
system of coupled oscillators on Z. We present this model immediately in the infinite volume setting. 

For each site n € Z, we associate an ordered pair (q n ,p n ) G R 2 and denote by X the set of all sequences 
x = {( < ZmPn)}nGZ- The system is said to be comprised of an infinite collection of oscillators, each situated 
at a site n £ Z with position g„gl and momentum p n G R. Here the state of the system is described by a 
sequence x = {(q n ,Pn)}neZ G X. 

The oscillators evolve in time according to the following coupled system of differential equations. For each 
n € Z and any t G R, 

(2.1) q n (t) = Pn (t) and p n (t) = e -(«»(*)-«»-i(*)) - e -(«»+i(*>-*»(*» , 

with initial condition given by some state of the system x 6 X. We note that the differential system (|2 . 1 [) 
corresponds to the following (formal) Hamiltonian H : X — >RU{oo} given by 

(2-2) fl"(x) = E Y + - <?n), 

ngZ 

where V(r) = e~ r + r — 1. For many sequences, the formal Hamiltonian may be infinite, however, it does 
formally generate the system of differential equations in (|2.1[) through Hamilton's equations. For classes of 
solutions with prescribed decay properties we refer to [42]. 

Existence and uniqueness of global solutions to (|2.1I) on certain subsets of X is well-known. Rather than 
address this directly, we begin by changing variables. Suppose that (|2.ip has a solution. For each n G Z and 
t G R, set 

(2.3) a n (t) = i e -C*H-i(*>-9»(t))/a a nd &„(*) = -ip„(t) , 

in terms of this given solution. This choice is commonly referred to as Flaschka variables as they were 
initially introduced in [12], see also [15] . Using (|2.1j) . it is clear that these new variables satisfy the following 
equations of motion 

(2.4) a n (t) = a n {t) (b n+1 (t) - b n {t)) and b n {t) = 2 (a 2 n (t) - a^t)) , 
and correspond to the following formal Hamiltonian 

(2.5) H (x) = J2 2b l + V (~ ln ( 4a «)) = E 2fe « + 4a « - 2 ln ( 2a «) - L 

nGZ nSZ 

In this work, we will concern ourselves mainly with the properties of the solutions of (|2.4j) . 

Consider the vector space M = £°°(Z, R) x £°°(Z, R) of pairs of bounded, real- valued sequences. We will 
write each x G M as 

(2.6) x = ({a„}„ eZ , {b n } n< zz) = {(a„, &„)}„ e z ■ 
M is a Banach space with respect to the norm 

(2.7) ||x|| M = max(||a||oo, Halloo), \\c\\<x> = sup \cn\ . 

n 

Given any initial condition x G M, global solvability of the system (|2.4|) is well-known (see e.g. [40], Theorem 
12.6). We will denote the Toda flow on M by <J> t , i.e., $t : M —> M is the function that associates an initial 
condition x G M to the solution of (|2.4[) at time t: $ f (x) = {(a„(i), 6„(t))} Tie z with $o( x ) = x - 

2.2. Estimating Toda solutions. More is known about the solutions of (|2.4j) on M. In fact, let us fix 
x e M and introduce operators L(x), P(x) : ^ 2 (Z) — > ^ 2 (Z) by setting 

(2.8) [L(x)/]„ = a n f n+1 + a„_i/ n _i + 6 n /„ 
and 

(2.9) [P(x)/]„ = a n f n+ i - a„_r/ n _i . 

For each x G M, the existence of global solutions imply that L($ t (x)) and P($t(x)) are well-defined for all 
tgl. Whenever the initial condition x is fixed, we will write L(t) — L($t(x)) and P(t) — P($t(x)) to spare 
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notation. Observe that for each x £ M, $t(x) £ M and so L(t) is a bounded self-adjoint operator with the 
operator norm satisfying 

(2.10) max(|| a (i)||oo,||6(i)||oo) < ||i(t)|| 2 < 2||a(t)|| 0o + ||6(t)|U. 

Similarly P(t) is a bounded skew-adjoint operator. Since P(t) is also differentiable with respect to t, it 
generates a two-parameter family of unitary propagators U (t, s) satisfying 

(2.11) — U(t, s) = P(t)U(t, s) with U(t, t) = t, 

for all pairs t,s £ R, see e.g. [10], Theorem 12.4 for more details. Moreover, a short calculation shows that 
P{t) and L(t) are a Lax pair associated to (|2.4p . i.e., 

(2.12) ±L(t) = [P(t),L(t)}, 
and therefore, 

L(i) = f7(i,s)i(s)?7(i,s)" 1 . 
From this fact, each x <E M satisfies the a-priori estimate 

(2.13) ||$t(x)|| M < ||L(t)|| 2 = ||i(0)|| a foralUeR. 

It is important to emphasize the fact that we work with real-valued solutions. In fact, while a local existence 
result can be established in the complex case, even periodic complex initial conditions can blow up in finite 
time (see e.g. [TJ] or [T5]). 

Our proof of the Lieb-Robinson bound for the Toda Lattice, see Theorem 12. 3[ makes crucial use of the 
following estimate describing the sensitivity of solutions to changes in the initial condition. 

Theorem 2.1. Let x £ M and /i > 0. There exists a number v = v(/i,x) for which given any n,m £ Z, the 
bound 

d 



oz 



OZ 



< g__ e -At(|n-m|-i;|t|) 

17 



(2.16) F n (i;x) 



(2.14) 

holds for all t £K and each z £ {a m , b m }. In fact, one may take 
(2.15) v= (l + VTtj ||L(0)|| 2 ( e " +1 + ^ 

Proof. Fix x € M. Global existence of solutions on M guarantees that for each x £ M and n £ Z, the 
function F n : R — > M 2 given by 

is well-defined. It is differentiable with respect to each z £ {a m , b m }, e.g. as a consequence of Lemma 4.1.9 
in pQ. When convenient, we will suppress the dependence of F n on x. Using the equations of motion, i.e. 
(|2.4|) . it is clear that 

The relation 

(2.18) j? n (t) = F n ( ) + ^ / D n>e {s)— F n+e (s)ds, 

Z Z I e| <1 2 

with 

(2.19) D n>e (s) 



(b n+1 (s) - b„(s)) S {e) a n (s)(-S (e) + <5i(e)) N 
4(a n (s)8o(e) - a n -i(s)S-x(e)) 

then follows immediately from (|2.1T|) . 

To complete our estimate, we introduce the following notation. For each v £ R 2 , we will denote by 

(2-20) v=(A and | v | = frj 
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Moreover, we will write 

(2.21) (^j < 0f\ if and only if \x\ < \u\ and \y\ < \v\. 

With this understanding, the uniform solution estimate (|2.13p . and (|2.19[) . it is clear that 

d — r\*\ 



(2.22) 



d_ 

dz 



F n (t) 



< 



|k(0) + ||L(0)|| 2 ]T (* D. 

|e|<l 



ds. 



where 
(2.23) 



n _ ( 2<5 (e) S (e) + 6 1 (e)\ 

1Je -[4(So(e) + 5- 1 {e)) ) ' 

Let us now consider the case that z = a m , i.e., 

d 

(2.24) 



In this case, iteration of (|2 .22[) yields 



^U(o) = (J) u») 



(2.25) 



d 



da, 



F n (t) 



OO 

* E- 

fc=0 



'W)Ut\) k 



k\ 



- E "' E S m+ei+ ... +ek (n)D ei ■ ■ -D ek ( Q j 

|ei|<l |e fc |<l V ' 



where we have set 
(2.26) 



* E 

fe= |n— m| 



(II^WIH*!)* 



D 



2 2 




D=Y,D e 

|e|<l 

Moreover, note that the remainder term in finite iterations of (|2.22j) converges to zero since -^F n (t) is 
continuous and thus bounded on compact time intervals. 

The eigenvalues of D are A± = 1 ± \/T7, and in terms of the eigenvectors v± given by 

(2.27) 



v± 



1 1 

u+ — V- 



it is clear that I ~ I = — i=v+ 

\ - / V / 2x/T7 2^17 

Taking the infinity norm, in K 2 , of both sides of (|2.25l) shows that 

Q 

(2.28) — F n (t) 

aa, 



< _ 
v 17 



_8_ ^ (X + \\L(0)\\ 2 \t\) h 

k— \n— m\ 



(2.29) 



A+||i(0)|| 2 . If c\t\ <\n- m\er^ +1 \ then by Stirling 

(c\t\) k ^ (C^l) |n ~ m| c e|t| < f C \A \ " e \n- m \ e c\t\ < g-^ln-ml-f |t|) _ 

' ^ Z^l IT) — in 1 1 \ j£ — jji I / 



Now, let /i > be fixed and set c = 

■ < 

fc! (n — m)! 

fc— |n— m| 

Otherwise < — fi(\n — m\ — ce M+ |t|), and so by (|2.28j) we have 
(2.30) '' 



da, 



■F n (t) 



In the case that z = b m , it is clear 



< 



e o\t\ < 



- /1 (| T ,- m |-c(e"+ I + i)|t|) 



that 



A_ 



A 



A short calculation leads to a prefactor of ^= in the analo 



We have proven the result 
Remarks 



'gue of (|2.28p , and this is less than the one above. 

□ 
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1. The number v in our estimate (12.14)) depends on the initial condition x only through the quantity ||L(0)||2- 
In fact, one can replace ||L(0)||2 with any uniform estimate on the solutions (uniform in n and t) as is 
clear in the bound from ([2TT51) to (|2"T2"2"|) . 

2. If \n — m\ > 0, the left hand side of (|2.14p is zero when t = 0. The estimate in (|2.14j) then shows that the 
solutions of the Toda Lattice at site n are insensitive (i.e. exponentially small) to changes in the initial 
condition at site m for times t satisfying v\t\ < \n — m|. For this reason, the number v is often called a 
bound on the velocity (i.e. rate) at which disturbances propagate through the Toda Lattice. Moreover, if 
\n — m| > 1, then the bound in (|2.28j) is at most linear in \t\, for small t. 

3. In general, the fact that solutions of the Toda Lattice have partial derivatives that decay exponentially 
in \n — m\ will be more important than any particular decay rate /i. For this reason, given an initial 
condition x G M, one can optimize the quantity v over all possible /i > 0. The graph of /(/x) = e M+1 + — 
is given below. 




(Mo,/(Mo)) 



(2.32) 



It is clear that the optimal /i is achieved when 







and moreover, this value is independent of the initial condition x. The unique solution of (|2.32[) can be 
expressed in terms of the Lambert W-function ([351 §4-13]) and is given by 

(2.33) /i = 2W(l/(2Ve)) w 0.47767 and f(fj, ) « 6.47622 . 

2.3. A Lieb Robinson bound for Toda. We can now formulate a Lieb-Robinson bound for the Toda 
Lattice. In general, an observable A is a function A : M — > C. Let us denote by A the set of all observables 
with well-defined first order partial derivatives. We will say that an observable A is pointwise bounded if 

(2.34) \\A\\ K = sup|A($ t (x))| < oo VxeAf. 

ten 

Take A^ 1 ' to be the set of those observables in A for which all of the first order partial derivatives are 
bounded and 



(2.35) 



u\\ 



E 





dA 




dA 




( 


da n 


+ 

X 


db n 


J 



< oo VxgM. 



An observable A e A is said to be supported on X C Z if the observables J^- and J^- are identically zero 
for all n € Z \ X, The support of an observable A is the minimal set on which A is supported, and we will 
denote this set by supp(A). Set Ao C A to be the set of all observables with compact support. For any 
t€M and A € A., the Toda dynamics, which we denote by at, is the observable- valued mapping given by 

(2.36) a t (A)=Ao$ t , 

where $ t is the Toda flow described at the end of Section |2~T1 above. In other words, for each x G M, tel, 
and A £ A, 

(2.37) [a t (A)](x)=A(* t (x)), 
which, since solutions are global on M, is a well-defined quantity. 

Example 2.2. The two most basic observables correspond to evaluation maps, i.e., to each n £ Z, we 
associate the functions A n and B n given by 

(2.38) An(x) = a n and B n (x) — b n for each xG M . 
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Each of these observables have support X — {n}. Moreover, the dynamically evolved observables 
(2.39) [a t (A n )](x) = a n (t) and [a t (B n )](x) = b n (t) 

correspond to observations over time. 

The modified Poisson bracket of two observables A and B is formally denned as the observable 

OA dB dA dB' 



(2.40) 



[{A,B}] (*) = jE' 



da n db n db n da n 



(x) 



where we have denote by -rf- 



>) 



gb — - -t^-. If, for example, either A or B has compact support, then the 
corresponding modified Poisson bracket is a well-defined observable. This quantity is of particular interest 
since it generates the Toda dynamics, i.e., for any A £ Ao 

d 



(2.41) 



dt 



a t (A) = a t ({A,H}) = {a t (A),H} , 



where H is the Hamiltonian (|2.5[) . Now, for any n,meZ,iel, and x € M, it is easy to see that 
(2.42) 



1 d 1 d 

{a t (A n ),B m } (x) = -a m -x- a n (t) - -a m - — a n (t) 

4 da m -i 4 da m 



using the basic observables defined in Example 12.21 In this simple case, Theorem 12.11 gives the bound 

2IIge.II 



\{a t {A n ),B m }(x)\ < 



e fj,(\n-m+l\-v\t\) _j_ e -fj,(\n-m\-v\t\) 



(2.43) 



< 



2 IHI°° n + e »\ e -»(\n-m\-v\t\) 



In general, we have the following Lieb-Robinson bound for the Toda Lattice. 

Theorem 2.3. Let x G M and \x > 0. There exist numbers C and v for which given any observables 
A, B € AS 1 ', the estimate 



(2.44) \{a t (A),B}(x)\<C NU ]T 





OA 




&4 




dB 




dB 




( 


da m 


+ 

X 


dbm 


X 


da n 


+ 

X 


db n 





,-/j,(\n-m\-v\t\) 



holds for all t G M. .Here C = -j^(l + eA w * s a s in Theorem \ 2. 1[ 
Proof. Substituting into (pT40|) . we find that 

d 



(2.45) 



{a t (A),B} ( X ) = i^a„. 



,,,95 d dB 
- — a t (A) ■ —a t (A) ■ - — 

da n n ' db n db n da n 



(x). 



Using the chain rule, it is clear that for any z G {a n , b n }, 



(2.46) 



a t (A) 



(*) = E 



at 



dA 

da m 



d_ 

dz 



att{A r , 



dA 
db m 



d_ 

d~z 



a t {B m ) 



(x) 



where we are using the notation from Example 12.21 Applying the triangle inequality and then Theorem 12. 1[ 
the bound 



\{a t (A),B}(x)\ < M^J2 



dA 



da r . 



a — a m {t) 
da n 



dA 



db n 



da n 



dB 



{ Moo ^ 



<9a n 



9L 



&4 



9 



&m(i) 



as 



< 



2||a|| 



(2.47) 



17 

2II0II 



E 



<9A 




<9A 




dB 


da m 


+ 

X 




J 


db n 



dh 

fj,(\n— m\ — v\t\) 



da n 



17 



"E 





dA 




<9A 




as 


( 


da m 


+ 

X 


db m 


J 


<9a n 



-/i,(|n— m| — -u|t| ) 
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readily follows. Note that the bound on e.g. 
triangle inequality proves (j2.44j) as claimed. 



follows from the argument in Theorem 12. II Another 

□ 



For many applications, the observables of interest will have disjoint supports. In this case, the bound in 
(|2.44[) can be stated as follows. For any x S M and /i > 0, there exist numbers C and v, as above, such that 
given any two disjoint subsets X and Y of Z the bound 

(2-48) \{a t (A),B}(x)\ < q|a||ooPI|i, x ||S||i, x e-' l(<i(x ' v ' ) - !j|t|) , 

holds for silt € M and any observables A,B<E A^ 1 ' with supports in X, Y, respectively. Here 

(2.49) d{X, Y) = inf{|x - y\ : x € X and y e Y} > , 



and (|2.48[) corresponds with the bound (|1.4p claimed in the introduction. Other norms on the observables 
under considerations are often useful. For example, consider 



(2.50) 



\dA\\* 







8A 




OA 




sup max 




da m 


7 


db m 








X 


X- 



for any A g A^ . Theorem 12.31 then shows that: For any x € M and fx > 0, there exist numbers C and v, as 
above, such that given any two disjoint subsets X and Y of Z, at least one of which being finite, the bound 



(2-51) 



\{a t (A),B}( X )\<4C\\a\\ 



|«9A|| x ||aB|| x e -MI^-"l-«l*D . 

m£X,neY 



holds for all t € M. and any observables A,B<E A^ with supports in X, Y, respectively. 

In any case, the bound in Theorem 12.31 demonstrates that each x £ M propagates no faster than some 
finite rate. We say that the Lieb-Robinson velocity corresponding to x g M is v(x, /io) where fig is the 
infimum defined in the remarks after Theorem 12. II 



2.4. On one-soliton solutions. Perhaps the most important collection of solutions to the Toda Lattice 
are the solitons, i.e. the solitary waves, see e.g. (40) [4TJ [43] . In fact, solitons can be considered as the stable 
part of any short-range initial condition since every such solution eventually splits into a number of stable 
solitons plus a decaying dispersive tail [22) . Since all involved quantities can be computed explicitly for the 
one-soliton solution we will use it as a test case and for our Lieb-Robinson bounds. 
Fix k > 0. A one-soliton solution of the Toda lattice is given by 

1 + exp [— 2nn ± 2 sinh(«;)t + S] 



(2.52) 



q n (t;±) = q-ln 



1 + exp [-2k(ti - 1) ± 2 smh(n)t + 5} , 

where q and 5 are real constants. Here q n (t;±) represents the position of the traveling wave. It describes 
a single bump traveling with speed ± sill M K ) anc j w idtli proportional to In other words, the smaller 

the soliton the faster it propagates. Changing 5 amounts to a shift of the solution and we will set 6 = 
for simplicity. From this explicit formula it is clear that, in contrast to harmonic models, the velocity of a 
solution to the Toda Lattice may indeed depend on the initial condition. 
In terms of the function 



(2.53) 

it is clear that 



(2.54) <?„(*; ±)=g- In 

and from Hamilton's equations, we also know that 
(2.55) 

In Flaschka's variables, we have that 
(2.56) an(*;±) = 



f±(x, t) = 1 + exp [— 2kx ± 2 sinh(«;)i] 
f±(n,t) 



/±(n-M) 



f .,s d f ± {n-U) f ± {n,t) 

p n (t;±) = —q n (t;±) = — — - . . 

at f±(n-l,t) f±{n,t) 



ly/f±(n-l,t)f±(n + l,t) __ J LU _^ 1 ( f'±{n,t) f ± (n-l,ty 



f±(n,t) 



and b n {t;±) 



2\f±(n,t) /±(n-l,t), 
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A short calculation shows that 

coshM , ,, , , , ,, , . , sinh(K) tanh(«) 
(2.57) supa„(0;±) =a (0;±) = and sup |6„(0; ±)| = |6o(0; ±)| = ^ ^. 

With x K denoting the initial condition corresponding to this one-soliton, we have proven that 

cosh(/c) /cosh(«) sinh(«) tanh(«) \ . „ , , . sinh(«) tanh(/c) 
(2-58) — ^ = max ^— ~^ —J < IWOIh < cosh(«) + U_ , 

using e.g. (I2.10[) . Since L(x K ) is self-adjoint its norm is equal to the spectral radius. Moreover, the spectrum 
is given by an absolutely continuous part [—1,1] plus the single eigenvalue ± cosIi(k) implying 

(2.59) ||Z(x«)||2 =cosh(/c). 

To see this last claim note that the one-soliton solution can be computed from the inverse scattering transform 
by choosing one eigenvalue A = ±cos1i(k) plus the corresponding norming constant 7 = (1 — e~ 2K )e s and 
zero reflection coefficient (cf. [43l Sect. 3.6] or [40l Sec. 13.4]) or by using the double commutation method 
to add one eigenvalue A with norming constant 7 to the trivial solution (cf. [40l Sect. 14.5]). 

As is clear from this calculation, the Lieb-Robinson velocity does provide a reasonable estimate on the 
actual velocity, at least for one-soliton solutions. 



3. Estimates on Perturbed Toda Systems 

In this section, we consider a class of perturbations of the Toda system for which locality results analogous 
to Theorem 12.31 still hold. We introduce these perturbations as follows. Let W : R — > [0, 00) satisfy 
W £ C 2 (R). Consider the formal Hamiltonian 

(3.1) H w = H + W n , 

nGZ 

where H is the Toda Hamiltonian, see (I2.5p . and for any initial condition x = {(a n , b n )} n ^z G ^0 = £°°(Z, R\ 
{0}) x £°°(Z, R), the observable W„(x) = W(ln(4a^)). Our choice of parametrization is motivated by the 
fact that the potential V for the unperturbed Toda Lattice, see (|2.2j) . is a function of q n +i — Sn = — m (4a 2 ). 
Since the choice a n = corresponds to taking the positions of the particles at sites n and n + 1 infinitely far 
apart, we will exclude it from our considerations. The formal Hamiltonian (|3.1[) corresponds to the following 
system of coupled differential equations 

(3.2) K(t) = <(t){K+i(t)-K(t)) 

bZ(t) = 2(<(t) 2 -<_ 1 (*) 2 )+iU*), 

where 

(3.3) R n {t) = \ [W/'(ln(4<(i) 2 )) - H/'(ln(4 a ^ 1 (i) 2 ))] . 

Local existence and uniqueness of solutions of (|3.2p corresponding to initial conditions x € Mq follows from 
standard results, [TJ Thm. 4.1.5]. Let us denote by the perturbed Toda flow, i.e., the function that 
associates initial conditions x 6 Mo with $7(x) = {(c^(t), (t))}, the solution of p. 21) at time t. Consider 
the set of initial conditions Mb = Mb(W) C Mo for which there exists numbers C±, C2 < 00 with 

(3.4) sup || $7 (x) || M < Ci and sup sup— < C 2 . 

teK tSK n£Z \a n (t)\ 

It will be shown in Appendix [A] that this M{, contains at least all initial conditions whose energy is finite 
under appropriate assumptions on W. 

For initial conditions x g Mb, we have two estimates on the corresponding Lieb-Robinson velocity. The 
first is obtained in Section 13.11 by simply arguing as we did in Theorem 12.11 The other, in Section 13.21 
achieves an estimate of the perturbed velocity (corresponding to x) in terms of the unperturbed velocity 
(corresponding to x) via interpolation. 
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3.1. Direct Bounds. Following closely the arguments in Theorem 12.11 we obtain the next result. 

Theorem 3.1. Fix W £ C 2 (R) with W" e /j > 0, and let x £ M b . There exist numbers C v 

C w (x, W) and v w = v w (fi, x) for which given any n,m £ X, the bound 



(3.5) 



d_ 

dz 



<(t) 



d_ 

dz 



Kit) 



< C ,w e ~' J (l n ~ m l -l;W l t l) 



holds for all t € M and each z £ {a m , b m }- In fact, one may take 
(3.6) 



v w = 1 



17 



±C 2 \\W"\ 



C x ( e^ +1 + - 



with Ci,C 2 from (I3T4I) . 

Proof. Fix x £ Mb. Global existence again guarantees that for each n £ Z, the function 

(3-7) FZ(t) = 

satisfies 

<(s) (b% + i(s)-b%(s)) 



(3.8) 



F:(t) = F:(o) + 



\2{aZ(s) 2 -a^) 2 )+Rn(s) 



ds. 



Since W is sufficiently smooth, the components of F™(t) are differentiable with respect to each z £ {a m , b m }, 
and the bound 

,|t| 



d 

—F w (t) 
dz " 1 j 



< 



|e|<l 



(3.9) 

follows as in f|2 . 22[) with 
(3.10) L>™ = 

Following the previous scheme, iteration (with z — a m ) yields 



DZ 



d 

dz 



els. 



2C 1 8 {e) C 1 (S (e)+S 1 (e)) 
(4C 1 + C 2 \\W"\\ oo ){5 (e) + 8- 1 (e)) 



(3.11) 



d 



da r , 



OO 

k— \n— m\ 



(2Ci|t|)* 



A:! 



(D w ) k 1 ' 



with 

(3.12) D 
Eigenvalues and eigenvectors of D w arc 
(3.13) 



a 



for a = 4 - 



C 2 \\W"\\ 
Ci 



1± Vl + 4a 
A± = and v± 



a 



and arguing as before, it is now clear that 

d 



(3.14) 
with 
(3.15) 

For z = b m , one similarly finds 
(3.16) 

This completes the proof. 



da r . 



< 



2a 



- fi(\n—m\—v w \ t\ ) 



VI + 4a 



2A+Ci e" +1 



db r 



< 



2X 



+ g-M(l«-m|-f w |t|) 



□ 
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3.2. Bounds Via Interpolation. The goal of this section is to prove a different bound on the Lieb- 
Robinson velocity corresponding to an initial condition x G Mb, see Theorem 13.31 below. The novel feature 
of this estimate is that it is explicit in the unperturbed Lieb-Robinson velocity of x. 

Before we state the main result of this section, we first indicate some further estimates on the unperturbed 
system which will be useful in proving Theorem 13.31 To start with, we prove another Lieb-Robinson type 
estimate for the Toda system, see Lemma 13.21 below. Afterwards, we introduce a quantity that is better 
suited for the iteration scheme which is at the heart of proving Theorem 13.31 see G M (fc) in (I3.37|) below. 
Lastly, we will state and prove Theorem 13.31 

We begin with the following lemma on second order derivatives of the unperturbed system. 

Lemma 3.2. Fix fi > and let x € M. There exists a number C = C{[i, x) > and a junction h, depending 
on fi and x, for which given any n,k,£ € Z, the estimate 



(3.17) 



max 



dzdbk 



~-a n {t) 



dzdb 



~-b n (t) 



< Ce-^-^e-^-^e^^hit) 



holds for all t € R. Here a n (t) and b n (t) are the solutions of (12.41) with initial condition x <E M, 



_d 

Dfc+i 96 



-gr— , z £ {ai, bi}, and the 



iber v 



is as in Theorem 



The function h grows at most exponentially. 



Proof. Fix x e M. As in Theorem 12.11 it is clear that the function F n (t) is well-defined for each n G Z and 
f € R. In addition, both (j2 . 18[) and (|2 . 1 9[) still hold with the choice z — bk- Taking a second derivative, we 
find that 



(3.18) 



& 



dzdbk 



bk \t\<iJ° ^ dbk ozdb k 



+e(s) ds 



for each z € {at, bi}, since 
The bound 

(3.19) 



dzdb k 



F n (0) = 0. 



dzdb 



< 



,-»\n-t 



17 



E 

M<i 



ob k 



ds 



^ r\t\ 
+ ||L(0)|| 2 E / D 

\e\<l J ° 



dzdbi. 



~-F n+e (s) 



ds 



with 
(3.20) 



D' 



(e" + l)<y (e) 6 (e) + S 1 (e] 
4(5 (e) + e^_ 1 (e)) 



and D e as in (|2.23j) . follows using Theorem 12.11 and the argument therein. A further application of Theo- 
rem [2TT] implies that 



(3.21) 



d 
db k 



< 



8 



17 



e t^ e ~fJ-\n+e-k\ e tivs 



(3.23) 



and therefore the first term on the right hand side of (|3.19l) can be estimated by 

(3.22) ?_ e V e -n\n-i\ e -v\n-k\ / e 2»vs ds 

Here we have used that 

e -li\n+e-k\ D' e ( l \ < e -^l"- fe l 
|e|<l ^ ' 

Let us introduce the notation 



2(e" + 1) 
4(e 2 " + 1) 



2(e^ + 1) 
4(e 2 " + r 



(3.24) 



C u = 



64 
17 f 



and y = 



2(e" + 1) 
4(e 2 ^ + 1) 
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We have proven that 
d 2 



(3.25) 



dzdb, 



r\t\ r\t\ 
< C /U e-^l"-^le-^l"-- fe l / e 2 ^ vs ds-y + \\L(0)\\ 2 / ^ D 
Jo Jo |e| < 1 



<9z<% 



ds 



Iteration now yields 
d 2 



dzdbh 



1*1 /•*! 



< 



C^||L(0)|P 2 / / 
J=0 Jo Jo 



(3.26) 



x ... e -p|«+eiH hej-^| e -/j|n+eiH he j — A: | ^-j^ jj^ 

|ei|<l |ejl<l 

and the convergence is guaranteed as before. 



D e 3 V 



A short calculation shows that 

-|t| 



(3.27) 



^vtj+t dt . +i ...^ 



i_ (2^|iir 



(2uw)j t± •« — ' m 

v ^ ' m=j + l 



whereas the bound e ^l z+e l < e^^'e M ' z l immediately implies that 

(3 28) e -M|n+ei+e 2 H hej-f| < e M|ei| . . . gM|ej | & -)J>\n-l\ 

This proves that 



(3.29) 

where 
(3.30) 



if 



dzdbk 



~-F n {t) 



< £lL e -iAn-t\ e -^\ 
~ 2/iv 



m— j + 1 



|e|<l 



2 e 2 ^ 1 + 1 

4(e 2 ^ + 1) 



(3.31) 



The eigenvalues of D are A± = 1 ± + 4(e 2 ^ + l) 2 and a convenient choice of eigenvectors are 

V± = (4^+ 1)) ' 
In terms of these, it is clear that 

2 (e M + 1) - A- , A+ - 2(e^ + 1) 

(3.32) y = yiv + + y 2 V- with yi = and y 2 = 



X. - A_ 



Inserting this into (|3.29|) . we find that 



dzdbk 



~-F n {t) 



< C » c -ix\n-i\ c -n\n-h\ 

~ 2[iv 



3=0 



A+ - A_ 

y\X 3 + v + + y 2 Aiu- 



m—j-\-l 



(3.33) 



j=o ^ 2fJ,v 



A-i 



E 

m=j-\-l 



(2nv\t\y 



It is easy to see that for any A, B e 

(B\t\r 



J2 AJ E 



(3.34) 



This proves Lemma 13.21 One may choose 



ml ' m! 

j=0 m— j+1 m— 1 



~ ( g | f |)m a . _ / ^ ( 6 CA-Dfl|*l - l) e *l*l A*l, 



B\t\e B W, 



A = l. 



(3.35) 



C=-f-Qy 1 \\\v + \\ 00 + \y 2 \\\v-\\ 00 ) 
2/j.v 



and by setting j3 = ^ L ^^ 2A+ — 1, we see that 



(3.36) 



h(t) 



(i^vm-i) if p ^o, 

\2^w|t| if /3 = . 
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Note that (3 depends only on /z, limt_>.o h(t) — 0, and moreover, if /3 < 0, then h(t) < |/3| 1 . For any /z, 
however, /i(i) < /ie r '*', where ft. = h(fi,x) and r = r(/i,x) > 0. □ 



Our proof of Theorem 13.31 will again use an iteration scheme. Due to interpolation, it will be necessary 
to sum certain terms over all integers. To this end, we make the following observation. For each /i > 0, the 
function 

P -Hfc| 

(3.37) G M (fc) 



satisfies the estimate 
(3.38) 



(l + \k\r 



/e: 



with 7 = 4£ feeZ (l + NT 2 - In fact, setting f(k) = (1 + \k\) 2 , it is clear that f{j + k)< 2(f(j) + f(k)), and 
thus 

Gfj,{j - ky'G^j - l)G^l - k) = f(j - k)f(j - 0" x /(i - fc)- 1 ^!^^-^'!-!'-^) 

(3.39) < 2(f(j-l)- 1 + f(l-k)- 1 ). 

The claim follows after summing over I. Clearly the choice of power 2 in the denominator of G M is merely 
for convenience; a factor (1 + |fc|) _1_<5 for any S > would suffice. Note, however, that there is no bound of 
the type in p. 381) for the exponential function with no inverse polynomial weight. 
As a final comment, it is clear that for any [i > and e > 0, the bound 

(3.40) e -(M+«0M < C B G„(|x|) 

holds for all x € K. The number C £ = sup^gjj (1 + |a;|) 2 e- E l x l < oo. 
We can now state the second result of this section. 

Theorem 3.3. Fix W € C 2 (R) with W',W" € L°°(R) ; £afce /i > 0, and Zet x e M 6 . for any e > 0, i/iere 
are positive numbers C — C(e), D — D(e, fi,x,W), and S = S(e, /i,x, W) such that 



(3.41) 



—b w (t) 



< CG„(\n - m\)e^ +E)vW \l + D (e m - 1 



holds for all t el and n, m G Z. //ere w = v(x, /i + e) is as in Theorem \2.1\ and z € {a m , o m }. 

In words, this result shows that for each /i > the perturbed Lieb-Robinson velocity v w (x, /x), corre- 
sponding to x s Mb, satisfies, for each e > 0, 

(3.42) " • ' ■ - > 



w w (x, /i) < 1 + - u(x, a* + e) + 



A* 



and, as is shown in the proof below, the dependence of 5 on W can be made explicit 

(3.43) S = 5 l \\W'\\ 00 +5 2 \\W"\\ 00 +8 3 
for some Si — 5j(e,/i, x) and i = 1,2,3. 

Proof. Fix W, x € Mfc, and tel. We begin by interpolating between the Toda and perturbed dynamics. 
As in Section |2~T1 denote by 

(3.44) a t (A)=Ao$ t and aJ(A) = Ao<S>J for any A £ A, 

the Toda and perturbed dynamics respectively. For any A € A, the observable-valued equation 

(3.45) aJ(A) - a t (A) = J ^-aJ(a t . s (A))ds 



is clear. Due to the form of the perturbation, 



-<(a t _ s (A)) = <({a t _ s (A),iJ w })-<(a t _ s ({A,^})) 



(3.46) 



j2a:(Wt- s (A),w k }) 
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and moreover, one finds that 

(3.47) {«t_.(A), W k } = -\w' k ■ J-at-^A) . 

1 ob k 

Combining these expressions, we have shown that 

(3.48) aJ(A) = a t (A) - ~ £ f < • J-a*-.^)) ds 

V db k J 

as a formal expansion. For specific choices of observables, the expression above can be estimated. 

Fix n £ Z and denote by A n and B n the observables introduced in Example l2.2l Consider the observable- 
valued 



(3.49) 



«r(A n ) 



and F n (t) = 



a t (A v 
a t {B n )J ' 



which for any given x £ M satisfies 



(3.50) 



[^T(t)](x) = 



= F%(t) and [F n (t)](x) = 



a n (t) 

bn{t) 



in terms of our previous notation. Our starting point is the equation 



(3.51) 



ajU-a t . s (B n ) i 



F n (t), 



ds . 



Since W is sufficiently smooth, the formula 



(3.52) 4-Jn(t) 



i)z 



dz 



where 
(3.53) 
and 
(3.54) 



D k (s;n) = 



^W^A^)-aj(-£- h a t - s (A n )) N 



D k ^(s;n) 




a t - s {A n )\ a 
a t - s (B n )) al 



w / d 



s \db e db k 



dbedbk 



a t - s (A n ) 
a t - s (B n ) 



readily follows. 

For any x £ M, the estimate 



(3.55) 



d ~ 

—F w (t) 
dz " 1 ' 



(x) 



< 



oz 
\\W'\\ 



(x) 

E 



E 



D k {s;n)^(s) 



I'l 



D k , t {s;n)^Ff(s) 
oz 



(x) 



(x) 



ds 



ds 



is clear. 

Using Theorem 12.11 the first term can be bounded by 

d 



(3.56) 



oz 



(x) 



Fn(t) 



< 



< 



-(/i+ e )(|n-m|-t)|t|) 



17 



where we have set u = u(x, /x + e) and used (|3.40|) . 
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For each x £ Mb, the matrix appearing in the second term satisfies 
(3.57) Wkten)]^ < C 2 \\W"\\ 



e l*+e e -(lJt+e)(\n-k\-v( s )Qt\- s )) 



17 

< C 2 \\W"\\ 00 -^=e' 1+e C e G f ,(\n - k\)e^ +e)v ^ tl ~ s) 
V17 



where we have denoted by 

(3.58) „00 = v($Y(x), fi + s) = (1 + VT7) ||L(*rW)lla + fa + ^) 

the unperturbed Lieb-Robinson corresponding to the initial condition $"(x). Since x £ 

(3-59) l|i(*rW)|| 2 < 2 |K(s)IL + ITOIL < 3d , 

and so the quantity in (|3.58l) can be estimated independent of s. Clearly, 



v = v(0) < swpv(s) := v* 



(3.60) 



The matrix in the third term can be dominated using Lemma 13.21 In fact, the bound 

(3.61) II MIL ^ 2Ce-^+ £ )l™- £ le-^ +e)|n - fc| e 2( ' i+e)l ' (s)(l *^ s) /i(|i| - s) 

follows immediately. Some comments are in order. First, the prefactor C from Lemma 13. 2[ appearing 
above, seems to depend on s through the velocity v(s); see (|3.35[) . If, however, one repeats the argument of 
Lemma RT21 and replaces the v(s) in (|3.27p with v* above, then the new prefactor is independent of s. Next, 
as discussed at the end of Lemma 13.21 the function h grows at most exponentially. In fact, in this particular 
application, it is clear that 

(3.62) h{\t\ - S) < ^ e 2(M+e)7^(|t|-^) 

with numbers h = h(x, fJ. + s) and 7 = 7(x, /i + e), each independent of s. This proves that 

(3.63) \\[D k As;n)\ (x)^ < CG,(\n - t\)G^\n - fc |) e 2(7+D(^K(|t|-«) . 
Putting everything together and suppressing the x-dependence, we have found that 



oz 



< B 1 GJ\n-m 



(3.64) 



e (M+=Ht| +fla y*G M (|n-*|) / e 

fcGZ J ° 

Bs E G^n-^G^n-k]) ( 



2(7+l)( A1 +eK(|t|-s) 



( M +e)«*(|t|- S ) 
0> 



ds 



< B 1 GJ\n-m\)e 



dz 

)(fj,+e)v*(\t\-s) 



ds 



ds 



Upon iteration, we find that 
d 



dz 



<B 1 G li {\n-m\)e a \ t \+B 1 Y, Bi E""E G ^ n ~ ' ' " G ^ ~ m|) 

j=l k x eZ kj&L 



1*1 rsi 
Jo 



where we have set 
(3.65) 

From (|3.38l) . it is clear that 
(3.66) 



e b(\t\- sl ) e b( sl - S2 ) . . . e b(.sj-i-Bi) e aB ida . ... dsi 

a = (ji + s)v and b = 2(7 + l)(/i + e)v* . 



E-E ^(In-AiD-.G^-mD^yG^ln-ml) 
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e KI*|-si) e K*l-*2) 



and the iterated integral can also be calculated: 
1*1 

(3.67) 
This shows that 
(3.68) 



Ksj-i-si) e ™ids r --d Sl 



,b\t\ 



I'l 



,(a-b) Sj d . . . dsi 



ob\t\ 



(a -by 



E 

k=j 



((a-b)\t\) k 



k\ 



<BiG^(|n-m|) e 



a\t\ , b\t\ 



Since v M + e < v*, it is clear that a — b < 0. In this case, 



3=1 



(3.69) 

follows from (|3.34[) and so 
(3.70) 



E 

fc=j 



fc! 



E 



Bj 



Bj 
a — b 



1~«a-b)\t\f 
k=j 



hi 



Bj 



<B 1 G Al (|n-m|)e a l t l ( 1 



Bj 



oBl\t\ _ Ja-b)\t\ 



^B 1 +b-a)\t\ 



B 1 



This proves (|3.41j) and completes the proof. 



□ 



3.3. Comments on these results. It is clear that results analogous to Theorem 13 . 1 1 and Theorem 13 .31 hold 
for more general, finite range potentials W. We stated the results as above for simplicity of presentation. 

It is also clear that analogues of Theorem 12.31 for more general observables, hold as direct corollaries of 
Theorem 13 . 1 1 and Theorem l3.3l Since the statements are clear, we do not rewrite them for the sake of brevity. 

4. Bounds for the Hierarchy 

In this section, we will demonstrate that the results from Section [2] also apply to the Toda hierarchy. We 
begin, in Section I4.1[ by introducing the Toda hierarchy using the recursive approach from [8] . We will use 
essentially the same notation as [401 and refer the interested reader to this text for further information. In 
Section l4~2l we state the Lieb-Robinson bound valid for the hierarchy, see Theorem 231 To prove this result, 
we first establish the crucial solution estimate, an analogue of Theorem 12.11 which shows that a solution of 
the hierarchy at site n has weak dependence on the initial condition at site m if \n — m\ » 1. The bound is 
explicit in terms of a quantity v r which, in particular, depends on the initial condition through the infinity 
norm of a specific matrix, see (|4.11j) . To make our bound even more concrete, we estimate the norm of this 
matrix in Lemma l4~3l below. Theorem 14.41 follows as in Section [2] 

We note that our results also immediately apply to the Kac-van Moerbeke hierarchy. This hierarchy can 
be viewed as a special case of the Toda hierarchy which is obtained by setting b n = in the even order Toda 
equations. As is discussed, e.g. in [27], this gives precisely the equations of the Kac-van Moerbeke hierarchy. 

4.1. The Toda Hierarchy. We introduce the Toda Hierarchy as follows. Fix r e N = N U {0}. Take 
cq = 1 and choose constants Cj £ K for 1 < j < r. For each x £ M, define sequences g(x,r) and h(x, r) 
componentwise by setting 



(4.1) 

and 
(4.2) 



g n (x, r) = J2 Cr-fgV+V (x) with # (x) = (S n , L(x)^„) , 



h n (x, r) 



3=0 



7 j C r-j 



^ +1 )(x) with ^')(x)=2a n <<5 n+l! L( X )^ n ), 



3=0 



with L(x) as in (|2.8[) . We now define a system of equations for the components of an unknown sequence 
x (t, r ) = {(a n {t,r),b n (t,r)}: 



(4.3) 



d n (t,r) = a n (t,r) (g n+1 (t,r) - g n (t,r)) 
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and 

(4.4) b n {t,r) = {K{t,r) - h n - x {t,r)) 

with initial condition x(0) = {(a„(0, r), b n (0, r)} = x£ M. As in the previous section, we have, for example, 
denoted by g n (t, r) = g n (x(t, r), r) to simplify notation. The system (|4.3[) and (|4.4|) is known to have global 
solutions (see e.g. [30], Theorem 12.6) for initial conditions in M. Varying r G No describes the Toda 



hierarchy. Let us denote by <& t r ^ the flow corresponding to the Toda hierarchy, i.e. ' : M — > M satisfies 
$( r) (x) = {(a n (t,r),b n (t,r)}, the solution of (|4~B1 and P~4l above, with $ r) (x) = x. 
The simplest example corresponds to r = 0. In this case, the system becomes 

(4.5) a n {t,0) = a n (t,0){b n+1 (t,0)-b n (t,0)) and 6„(i, 0) = 2{a n (t, 0) 2 - a„_i(t, 0) 2 ) , 

which is, of course, the Toda system (|2.4I) . 

4.2. The Lieb Robinson Bound. The above choice of sequences g and h guarantee that the Lax-formalism 
of the Toda Lattice, see (|2.8[) - (|2.12j) . still holds for the hierarchy. In fact, set 



(r) 



(4.6) P(x,r) = Y J Cr-jP {j+1) (x) with P U) (x) = [L(xy>]+ - [L(x) J ]_, 

i=o 

where [A]± denote the upper and lower triangular parts of an operator with respect to the standard basis 
drain) = <5m,n (with 5 m>n the usual Kronecker delta). It is known, see e.g. [40 , that the Toda hierarchy is 
equivalent to the Lax equation 

(4.7) 



±L(t) = {P(t,r),L(t)}, 



where we have set L(t) — L(${ (x)) and P(t,r) = P($[ r ' ) (x), r) again to ease notation. It is easy to see 
that the operator P(t, r), which is of order 2r + 2, is skew-adjoint and differentiable. Like before then, there 
exists a unique unitary propagator U^ r '{t, s) for P(t,r). It follows from the Lax equation that 

(4.8) L(t) = U^ r) {t, s)L(s)U {r '>(t, s)' 1 

implying again an a-priori estimate 

(4-9) maxdlafttOIL ,||6(*,r)||J < ||L(f)|| a = ||L(0)|| 2 . 

For later use we record the following structure. 

Lemma 4.1 ([23 ). For each x G M and any integer j > 1, the sequences g^\x) and h^(x), as defined 
in (|4.ip and (|4.2[) . have components that are homogeneous. In fact, they are sums of monomials of the 
components of x with degree j and j + 1, respectively, which have the form: 

fk-i \ 

b n+k + R{n + k — l,n — k + 1)+ 



+ 



n a 2 n+i 

1=0 

flat 

fc-2 

II a 2 n +i 

1=0 



bn—k 



n+fc-1 

+R(n + k - 2, n - k - 



k— 1 

2 £ b n -> 

1=0 



n 



3 = 2fc -f 

fc-2 

2&n+fc-l b 



and 



j#>(x) = 



fc-1 



2 n 



•=0 



'n+fc 



6 2 

"n+fc 



fc-1 V 
2&n+fc £ &n+i) + 

fc 

R(n + k - 1, n - k + 1) + 2 Jl j = 2fc + 1, 



fc-i 



2 n 



n+£ 



=0 
fc-1 

2 n <- 

1=0 



R[n- 



1=0 

k — 1, n 



k + 2)- 



fc-2 

b n +i + b n -k+i +2J2 b n - 

1=0 



3 = 2k, 
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for j > 1. Here R(n, m) denotes terms which involve only at and bg with m < £ < n and we set R(n, m) = 
if n < m. 

The result below is an analogue of Theorem 12.11 for the Toda hierarchy. 
Theorem 4.2. Fix x € M, fi > and r £ Nq. For any n, m E Z, the bound 



(4.10) 



d_ 

dz 



a„(t,r) 



0_ 

dz 



b n (t,r) 



< e 



-/i(r|n-m|/(LSJ+l)l-iv|t|) 



holds for all i 6 1, where z G {a m , b m }, 

(4.11) v r = v r {x,fx) = |p(r)|U|L(0)|| 2 + j^j , 

and D(r) depends on x ; r, and the numbers ci, . . . , c r . 

Proof. Without loss of generality assume t > 0. For each n € Z, define the function F n : R — > R 2 by 

/ o„(i,r) 



(4.12) 



Fn(*,r) 



Since r will be fixed for the remainder of the argument, we will drop it from our notation. Using the equations 
of motion (|4.3p and (I4.4[) . it is clear that 

fi 



(4.13) 



F n (t) = F n {0) + 



a n( s ) (g n +l(s) - 9n{s)) 
h n (s) - /lft_l(s) 



ds 



Observe that by Lemma |4~T1 for each < j < r, the quantities a n (s)(g^P (s) — 9n (s)) and hn +1> (s) ~ 

hn-l\s) are homogeneous polynomials (of degree j + 2) in the variables a n+e {s) and b n+e (s) for |e| < |_§J +1- 
Differentiating with respect to z € {a m , & m } we get 

<9z 



(4.14) 



F„(t) = ^F„(0) + £ Cr _, Yl J* D^ e (s)^-F n+e (s)ds, 



3=0 |e|<L4j+l 



where the entries of Dn\(s) are homogeneous polynomials of degree j + 1. In particular, the following 
estimate, analogous to the bound in (|2.22[) . holds 



(4.15) 



\D^ e (s)\<\\L(0)\\i +1 D^, 



an explicit formula for appears in (I4.27[) below. As a result, it is clear that 



(4.16) 



0_ 

dz 



F n (t) 



< 



F n (0) 



Ei^-m^nf 1 E f D P 



j=0 \e\<[H+l 

Upon iteration of the above inequality (in the case that z = a m ), we find that 

„ oo j, r r 



d_ 

dz 



F n+e (s) 



ds . 



da, 



-F n (t) 



(4.17) 



< 



< 



fe! 

k=o ji=o Jk =o 

x E "' E ^m+ei + -+e«, 
|ei|<L^J+l M<L^J + 1 

(\\L(0)ht) k 

3i=0 



E 

fe=nn-m|/(L5J+l)l 



E 

fc=n»»-fn|/(L5J+l)l 



(\\L(0)\\ 2 t) k k (l 



kl 



E icr-jsmM 1 - E \cr- jk \\\L(o)\\i"D^...D^ (I 

h=o 
D(r) 



jk=0 
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where we have set 



(4.18) = J2 D e ] and D ( T ) = ^2\ C r-i\\\ L (°)\\i D(i) 

\e\<li\+l 3=° 



A similar estimate holds for the case z = b m . 

Taking the infinity norm one obtains that for any \x > 

(»«„nM, lk 



(4.19) 



az 



where we have set v r = WD^)]^ ||£(0)|| 2 (e M+1 + ^) 



fc=r|n-m|/(L5J+l)l 

< e -A»(rh-H/(L5J+i)l-«r|*|) 

□ 



We now provide a rough, but explicit, estimate on the velocity corresponding to the Toda hierarchy. 

Lemma 4.3. Fix x G M , fi > and r € No- TVie velocity corresponding to the Toda hierarchy, see (|4.10j) 
and (|4.11j) satisfies 

(4.20) w r < 8 + -) E |cr-j|||£(0)||^ +1 C7 + 2)3 J 

v ^ 7 j=0 

Proof. It is clear from (14.111) that we need only estimate the quantity HZ^r)!^, with D(r) as defined in 
(|4.18[) . To see this, first recall that for each 1 < j < r, the quantities a n (s)(g^^ (s) — 9n (s)) an d 
hn +1 \s) — h^P(s) are homogeneous polynomials (of degree j + 2) in the variables a n+e (s) and b n+e (s) for 
| e j < [§J + !• We need an estimate on the number of terms in each of these polynomials. 

Both sequences g^ +1 ^> and h^ +1 \ see (|4.1[) and (|4.2[) . have components defined in terms of the operator 
L(x) given by ()2.8|) . In terms of the shifts S ± on £ 2 (Z), i.e. (S ± f) n — fn±i, one can write 

(4.21) L(x) = aS + + a~S~ +b. 

Here a and 6 are regarded a multiplication operators and a 1 * 1 is the multiplication operator given by (a ± /)„ = 
a n ±ifn- An upper bound on the desired number of terms can be obtained by taking both a and b to be 
constant sequences, i.e., calculating 

<«*> -<«.,(*♦+ 5- +iT«.)- £ ( J DG/2 

fc— 0, k even 

and 

(4.23) ^-(W + ^ + iT*.)- £ ( J ' 1 ( (fc + \ V2 ) 

fc=0, ft ocid v 7 vv ;/ 7 

The bounds t/ j+1 ) < and < 3 3 readily follow. 

Next, using this estimate, we can expand 

2^0 + 1) j + 2 

(4-24) a n {8)(gV+t\a)-gH +1 \,j)= £ II 

fc=i ^=i 

and 

(4-25) ^• +1 )( S )-^ 1 1) ( S ) = 2 E II<W M 

fc=i f=i 

where d„+ e £ {±a„ +e (s), ±6„ +e (s)}. Here we have inserted a 2 because of the definition of h^ +1 \ see (|4.2D . 
This form enables us to determine the entries of the matrix Dn} e (s). In fact, differentiation of the left hand 
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side of (I4.24p , with respect to z £ {a m , b m }, yields 
(4.26) 



2^W+!) j+2 j+2 2r,< 3 + 1) j+2 j+2 

fe=l o=l | e |<|_2j + i fc=l o=l 



and an analogous formula holds for the partial derivative of the left hand side of (|4.25l) . These expressions 
determine the matrix entries of Dn}e{s). The bound in (|4. 1 5[) now follows with 



(4.27) 

Summing on e yields, 
(4.28) 



|e|<LiJ+l 

where we have used r)V +1 > < 2£W +1 ). This shows that 



= 8(j + 2)C°' +1) 



(4.29) 



and we are done. 



P(r)|U < 8^ |c r _ 3 |||L(0)|| J 2 (j + 2)3* 



□ 



We end this section with an analogue of Theorem 12.31 for the hierarchy. Let otf ^ denote the dynamics 



corresponding to the Toda hierarchy, i.e., a { t r) (A) = Ao & t r) for all A £ A. 

Theorem 4.4. Let r £ No, x € M , and /J > 0. There exist numbers C r and v r for which given any 
observables A,B£ A^ x \ the estimate 



(4.30) \{a\ r \A),B}(x)\<C £ 





OA 




BA 




8B 




dB 




( 


da m 


+ 

X 


3b m 


J( 




+ 

X 


db n 


J 



,-|*(r|n-m|/(L5J+l)l-Ur|t|) 



holds for all isK. TJere 4C r = ||a||oo(l + e M ^^ L S J + 1 ) ^) Vr j s as jj-j Theorem 

Proof. The estimate (|4.30j) follows as in the proof of Theorem 12.31 using the results of Theorem 14.21 as 
input. □ 

5. Results for the Perturbed Hierarchy 

The purpose of this section is to demonstrate that the methods from Section [3] also apply to the Toda 
hierarchy. Many of the proofs follow closely the previous arguments, and so we only sketch the details. 

To introduce the relevant class of perturbations, we first recall the Hamiltonian formulation of the hier- 
archy, see e.g. Section 1.7 of [15]. For each r £ No, consider the following formal Hamiltonian 



(5.1) 



r + 2 



t^ij^o I ' rodd ' 



with g ( k 3+2) as defined in (|47Tj) . 



(5.2 



It can be shown that H r generates the Toda hierarchy in the sense that for any A £ Ao, 

±c£\A) = c£> ({A,H r }) = {a[ r \A),H r } 



where a t is the dynamics associated with the Toda hierarchy as introduced before Theorem 14.41 In fact, it 
is shown in e.g. [HI Thm. 1.71] that 

dH r 



(5.3) 



OH 

-(x) = 4g n (x,r) and A n — — -(x) = 4/i n (x,r) 

UQrr OO n 



and therefore the evolution equations 



(5.4) 



(A ri ) = a t (r) ({A„,ff r }) = 



dt 



8H r 

Bb n+ 



dH r 
db n 
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and similarly, 
(5.5) 



follow; compare with (|4.3|) and (|4,4|l 

Note that for r = 



dHr 



dH r 

On-1 



(5.6) 



ff Q (x)=£(2&2+4a 2 fc -l) 



is different from (|2.5|) . However, the equations of motion (12.41) are the same as above. 

We can now introduce perturbations as in Section [3j Fix r £ N and let W 
W £ C 2 (M). Consider the formal Hamiltonian 



[0, oo) satisfy 



(5.7) 



where H r is as in (|5.1|) above and W n is the observable with W n (x) = VF(ln(4a 2 )). 
The equations of motion corresponding to H™ are 

(5.8) <(t, r) = <(t, r) ( 5n+1 (i, r) - g n (t, r)) 
and 

(5.9) bZ(t,r) = fc„(t,r) - h n ^(t,r) + R n (t) 
where 



(5.10) 



= 5 [TT(ln(4<(i,r) 2 )) - W> (ln(4aj_ 1 (t, r) 2 ))] 



Again, local existence and uniqueness of solutions of (|5.8[) and (|5.9p . corresponding to initial conditions 
x € Mo, follows from standard results, [TJ Thm. 4.1.5]. As before, let us denote by <J>™ r the perturbed flow of 
the Toda hierarchy, i.e., the function $^ r (x) = {(a™(i, r), r))}. Our arguments apply to the set of initial 
conditions Mf, ir = Mb <r {W) with bounded trajectories, i.e., for which there exists numbers C±, C2 < 00 with 

(5.11) sup \\$Y r {x)\\ M < Ci and supsup — \ < C 2 ■ 

For initial conditions x 6 Mb, r , there are results similar to the two main estimates from Section [3] We 
first state an analogue of Theorem 13. II 

Theorem 5.1. Fix r £ No, W £ C 2 (tt) with W" € i°°(R), and x £ M 6>r . For each fj, > 0, there exist a 
number v™ = u^(/i,x) for which given any n, m £ Z, i/ie estimate 



(5.12) 



TOO 



< e 



- M (nn-m|/(LfJ+l)l-<|4|) 



ZioZds /or i £ R. _ffere z £ {a m , b m }. 

Proof. We will follow closely the proof of Theorem 13.11 using Theorem 14.21 as input. Take t > and, as 
before, introduce 

(5-13) ^ P) = (wr))- 

We will suppress the dependence on r. Since W is sufficiently smooth, F™ is differentiable with respect to 
z £ {a m , 6 m } and the bound 



(5.14) 



oz 



< 



r+1 „t 



"a - -^ra+e( s ) 



j=0 | e |</3 3 ' 

follows as in Theorem 14.21 Here we have denoted by 

(5.15) ctj = \cr_j\Ci +1 for < j < r with a r+1 = C 2 \\W"\ 
and 

(5.16) fa = [j/2\ +1 for < j < r with (3 r+1 = 1 . 



ds . 
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Moreover, for < j < r, we have taken as in the proof of Theorem 14.21 and set 



(5.17) 





1 



for e = —1, and D\ 



(r+l) _ fO 







Taking z — a m and iterating yields 

d 



(5.18) 

with 
(5.19) 



dz 



< 



E 



fc=r|n-m|/(L5J+l)l 



r+l 



= J2 a 3 DU) and DC,,) = E D e j) 
3=0 |e|<ft 

A similar estimate holds for the case z = &„,.. 

As in the proof of Theorem 14.21 the bound (|5.12[) now follows with 

1' 



(5.20) 

Note also that 
(5.21) 



\Dl 



Halloo < ICr-il^P^IU + ^H^'ll 

3=0 

and so the estimate from Lemma 14.31 applies here as well. 



□ 

The interpolation argument proven in Section [3] generalizes to the hierarchy as well. To see this, we first 
prove an analogue of Lemma [ 



Lemma 5.2. Fix r G No, (J, > 0, and let x G Mq. There exists a number C = C(r, fx, x) > and a function 
h, depending on r, /i, and x, for which given any n,k,£ G Z, t/ie estimate 



(5.22) 



dzdbk 



~-a n (t,r) 



dzdb 



~-h n (t,r) 



< C*e^ I ( n " _£|/(L -5 J+1)1 )e _Al ( nn_fe|/(L 5 J+1)1 )e 2AI1, '' |t| /i(t) 



/io/cfc /or aZ/ i G 1. Here a n {t,r) and b n (t,r) are the solutions of (|2.4j) wzt/i initial condition x G M, 
- -M-, z G {a£,&£}, and the number v r is as in Theorem \2.1\ The function h grows at most 



db k db k + 1 

exponentially. 

Proof. Again, since this proof follows closely the arguments of Lemma 13.21 we will only sketch the details. 
Keeping with the previous notation, it is clear that 

d 2 

(5.23) 



dzdb k \e\<W2\+l Ja ^ V dh J 



3=U |e|<Lj/2j+l 

Differentiation with respect to z produces two terms. We begin by estimating the first term, i.e., the one 



that contains -§-Dn}e{s). Using Theorem 14. 2 [ it is clear that if |e| < \J/2\ + 1, then 



d 



(5.24) 

For convenience, let us denote by 
(5.25) 

Using (|5.24j) . one can show that 

8 



db k 



< e^e 



m -M(n«-fei/(LiJ+i)i-^| S |) 



(5.26) 

where 
(5.27) 



E 

M<|j/2J+1 



I'l 



d 



E(x) = e-^Cn^l/CLfJ+i)!) . 

jj(fJ,)E(n - k)E(n jf e 2 ^ s ds Q 
7i ( M )=8e 2 lL(0)||i(j + l)(j + 2). 



ds < 
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This proves that 
d 2 



dzdbk 



r-F n {t) 



< 



Cp(r)E(n - k)E(n - 1) J e 2 ^" s ds Q 



(5.28) 

where 
(5.29) 



|L(0)|| 2 ^| Cl ._,|||L 



E 

|e|<U/2J+l 



1*1 



dzdb 



ds 



3=0 



and De is as in (|4.15p . Iteration, as in the proof of Lemma [3~2l yields 



(5.30) 

where 
(5.31) 



i:) 2 



dzdbk 



F n (t) 



< ^^±l E (n - k)E(n -£)^T 
2[iv r ^ 

r m— 



\mh e 2, 

2fj,v r 



g ^r\t\Y b{rf 



p—m-\-l 



D(r)=^|c r _ i |||L(0)g J2 D e ] 

J=0 |e|<Lj/2j+l 

Taking oo-norms, and using (|3.34[) , we see that 

d 2 



(5.32) 
where 
(5.33) 



dzdbk 



^F n {t) 



C (r) 

C = and with p = 

2fj,v r 



< CE(n - k)E(n ~ £)e 2tiv ' w h(t) 



2/i« r 



h is as in (I3.36[) . This proves (|5.22[) and we are done. 



□ 



Given Lemma I5T21 above, the analogue of Theorem 13.31 follows with only minor modifications. We state it 
below. 

Theorem 5.3. Fix r eN ,W E C 2 (R) with W, W" € L°°(M.), and let x € M b>r . For each ^ > and any 
e > 0, there are positive numbers C = C(r, e), D = D(r, e, fj,, x, W), and 6 = 8{r,e, fj,,x,W) such that 



(5.34) 



< CG^Qn - m\)e^ +e >r\t\ [ x + D Lfi\t\ _ x 



holds for all tel. Here v r — v r (fi + e,x) is as in Theorem \4.2\ [i r — /V(L§J + 1); an d z G { a mi b m }. 
Proof. Interpolating as before, it is clear that for any x E M, the bound 



(5.35) 



d - 

—F w (t) 

dz " ( ) 



(x) 



< 



dz 

. \\W\L 



(x; 



E 



E 

fcez 



D k (-s;n)^-F^(s) 



\t\ 



D M (s;n)^J7(s) 



(x) 



(x) 



ds 



ds 



follows as in the proof of (|3 . 55[) : here all quantities depend now also on r, but we have suppressed this in 
our notation. 

With Theorem 14.21 it is clear that 



(5.36) 



(x) 



^Fn{t) 

dz 



< e -O+<0(rin-m|/(L§J+l)"|-T>r|t|) 

< C £ , r G^(\n-m\)e^ +e)vM , 
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where we have set 

(5.37) C £ , r = sup(l+a;) 2 e- e (^/(L§J+i)l) ; fi r = and v r = u r (x, fi + e) . 

x>0 L2J + 1 

Similarly, for each x G r , the matrix appearing in the second term satisfies 

(5.38) II [£*(*;«)] MIL < C , 2 ||W / || 00 e^ +s e- ( ^ +e) (ri n -*l/(L5J+ 1 )l-^W(l*l-»)) 

< C 2 }\W"\\ 00 ei* +s C e , r G IM .{\n - k\)e^<^-^ 

with 

(5.39) v r (s) = v r ($7r( x )i A* + e ) an d v*—aupv r (s). 
Lastly, the bound 

(5.40) llpMfan)](*)IL <CG Mr (|n-£|)G Mr (|n-fc|)e 2 ^+ 1 )^+ e K^I- s ) 

follows as in the proof of (|3.63[) ; with (possibly) different values of C and 7. 

Iteration yields (|5.34p as in the proof of Theorem 13.31 □ 

6. Locality bounds for more general initial conditions 

In this section, we return to the class of perturbations considered in Section [3l see also Section [5] The 
goal here is to prove a locality estimate for more general initial conditions. Let us recall the basic set-up. 
Fix W : K -> [0, 00) satisfying W G C 2 (R) with W', W" G L°°(M). The formal Hamiltonian is given by 

(6.1) H w = H + J2 W n , 

where H is the Toda Hamiltonian as in (I2.5[) . and for each n, the perturbation W n is taken to be the 
observable W„(x) = IF(ln(4a 2 )). The corresponding equations of motion are 

(6-2) aZ(t) = aZ(t)(bZ +1 (t)-bZ(t)) 

bZ(t) = 2(al{tf^al_ l {tf)+R n {t) 

where 

(6.3) R n (t) = \ [VF'(ln(4<(t) 2 )) ~ ^'(ln(4a^ 1 (t) 2 ))] . 



As we discussed before, local existence and uniqueness of solutions of ()6.2[) corresponding to initial conditions 
x G M follows again from Theorem 4.1.5 in pQ. For x G Mo, let us again denote by <J>r( x ) = {( a n(*)i K (*))} 
the perturbed Toda flow. If we set, L{t) = L($™(x)) and P(t) = P($™(x)) with L and P as in and 
(|2.9[) . it is easy to check that 



(6.4) |i(t) = [P(t),L(t)]+JZ(i), 
where the multiplication operator R(t) : £ 2 (Z) —> is defined by 

(6.5) [R(t)f] n =Rn(t)f n . 
From this definition, it is clear that 

(6-6) ||JZ(f)|| 2 < HW'lloo, 

and hence the operator norm of R(t) is bounded uniformly in both t and the initial condition x G Mo. 

One can check that Pit) still corresponds to a family of unitary propagators, which we denote by U(t, s). 
A short calculation shows then that Lit) = U(t, s)* L(t)U(t, s) satisfies 

(6.7) j t L(t) = Uit,s)*R(t)U(t,s). 

In this case, 

(6-8) ||i(t)|| 2 - ||L(i)|| a < ||L(0)|| 2 + / * ||fl(*)|| a d a < ||i(0)|| 2 + ||W"|U|*| . 
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Since the bound 

(6.9) max(sup|<(0|,sup|TO|) < \\L(t)\\ 2 

n n 

holds, (|6 .8[) produces a linear bound on the growth of solutions. From this, the existence of global solutions 
follows from Proposition 4.1.22 in [T]. 

A result analogous to Theorem 12.11 follows . 

Theorem 6.1. Let x = {(a n ,b n )}nez £ Mq with a* = inf„ |a„| > 0. Then for any fi > and n,m £ Z the 

estimate 



(6.10) 



max 



d_ 

dz 



d_ 

dz 



Kit) 



< max 1 



- fi(\n—m\—v w (t)) 



holds for all t £ R. Here z £ {a m , b m } and v w (t) is an explicit cubic polynomial with i> w (0) = 0. 



Proof. Fix x £ Mo with a* > 0. Global existence on M guarantees that for each n £ Z, the function 
F™ : K -> R 2 given by 

dn«(t) 2 ) N 
*C(<) 

is well-defined for all tel. In fact, since a* > 0, the equations of motion ensure that 

rt 



(6.11) 



TO = 



<(*) = < (0)exp (6- x {s) - bZ(s))d, 



(6.12) 



preventing a singularity in the logarithm. It is then clear that (|6.2p implies 



V<(s) 2 -<-l( S ) 2 + |iin(s) 



ds . 



(6.13) F„ w (0=F„ w (0)+2 
For any z £ {a m , &,„}, a relation similar to (j2.18[) holds, i.e. 

(6-14) ^i^(t) = |^(0) + 2 E / j^„ w + e(s)^, 

where 



|e|<l' 







(6.15) 



The bounds in MM and (I6T91) show that 



w 



Si(e) - S (e) s 
5-x(e) 



(6.16) 

and so we find that 



<M 2 <(iiL(o)ih + n»"iw»ir , 



< 



% r\t\ 
E / 9(s)Dl 

\e\<l J ° 



(6.17) 

holds for any (£R, Here we have denoted by 

■> WWW 

(6.18) 5 ( s ) = i + (||i (0 )|| 2 + ||^'|| oo | s |) 2 + ^^ and D v 



<9 

9i F ™ W+e(s) 







5 (e) + 5_i(e) 



*i(e) +S (e) 




Iteration (with z = a m ) yields 
(6.19) 



(6.20) 
(6.21) 



* E E ••• E 

fe=|n-m| |e x |<l | e fc | < 1 

9 



<9z 



n+eiH hefc 



(0) 



dt]~ ■ ■ ■ dt\ 



< 



E 



h(t) k 



i°"»i*=~h fc! 



(D w )' 



0/ ' 
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where hit) = 2 jj 4 ' g(s)ds and 

(6.22) £> w = \ £ DJ = (I 

In the case that z = b m , it is clear that 
(6.23) 



|e|<l 



1 




d 

—F w (t) 



< t fW(: 



k=\n—m\ 



Taking the infinity norm, we get for any fi > 

d 



3.24) 



where v w (t) = h(t) 



d. 



■F w (t) 



< max 1 



E 



k—\n—m\ 



< max [ 1,— ] e-Mdn-mi-^ct))^ 
a* 



□ 



Again, everything extends to the hierarchy. Since the statements are clear, we do not rewrite them for 
the sake of brevity. 

Acknowledgements: The authors would like to thank Bruno Nachtergaele for suggesting the problem 
discussed in this article and for several productive discussions. Also, the authors would like to acknowledge 
the hospitality of the Erwin Schrodinger Institute in Vienna, Austria where parts of this paper were discussed. 



Appendix A. Existence of Bounded Solutions for Hamiltonian Systems 
In this appendix we want to look at a general Hamiltonian system with nearest neighbor interaction: 

(A.l) H(x) = E ( f + V (^+i In)) , 

where x = {{p n , q n )} n& . Let V € C 2 (R) with V{x) > and V(0) = V'(Q) = such that is a fixed point 
of the system. Since we want to obtain bounded solutions we will assume that our interaction potential is 
confining in the sense that V(x) +oo as x — >• ±oo. Since the uniform motion q n (t) = q a +p°t (with q°,p° 
some real constants) of the system is unbounded we switch to relative coordinates r n = q n +i ~ q n in which 
the equation of motions read 

(A. 2) r n =p n+ i-p n , p n = V(r n ) - V'(r n -i)- 

We will consider this system in the Hilbert space X = l?(TL) x £ 2 (Z). 

Theorem A.l. Suppose V G C 2 (R) such that is a unique global minimum with V(0) = V'(0) = 0, 
V"(0) > and V(x) +oo as \x\ -> oo. 

Then the system (|A.2[) has a unique global solution in X for which the energy 

(A.3) H{pjr) = J2(f + V(r n )\ 

is finite and conserved. This solution is C 1 with respect to the initial condition. Moreover, is a stable 
fixed point and all solutions satisfy || (p(i), r(t))|| 2 < C as well as || (p(<), r(t))|| oc < C, where the constant C 
depends only on the initial condition. 

Proof. First of all note that by our assumption on V we can find constants c# and Cr for every R > such 
that \V(x)\ < C R x 2 , \V'(x)\ < C R \x\ and \V{x)\ > c R x 2 for \x\ < R. 

In particular, for r G £ 2 (Z) with ||r||2 < R we have ||V'(r)||2 < Cfi||?"||2 and it follows that the map 
r n t— > V'(r n ) is C 1 on £ 2 (Z). Since the same is true for the shift operator x n i— > x n —%, our vector field is 
C 1 and local existence and uniqueness follow from standard results [1] Thm. 4.1.5]. This also implies that 
the flow is C 1 with respect to the initial condition [TJ Lem. 4.1.9]. Moreover, \H(p,r)\ < Cij||(p,r)||2 implies 
that H is finite on X and a straightforward calculation shows that it is conserved by the flow. 
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Moreover, H(p, r) > c\e for r)||2 = £ and e < 1. Hence [TJ Thm. 4.3.11] shows that is a stable fixed 
point. 

Finally, if V[x) — > +oo there is a constant Me such that |x| < Me whenever |V(af)| < E. Hence setting 
E = #(p(0),r(Q)) we have H(p(t),r(t)) = E implying \\p(t)\\ 2 < V2E and IK*)^ < M E . But this implies 
\\r{t)\\ 2 < cllJ\V(r)\\ 2 < \[E. Hence our vector field remains bounded along integral curves on finite t 
intervals and hence all solutions are global in time by [TJ Prop. 4.1.22]. □ 

Note that, using q n (t) = q n (0) + L p n (s)ds, for our original variables we get 

(A.4) ||g(i)||a < ||<z(0)|| 2 + C|t|, ||g(t)|U < ||g(0)||oo + C\t\. 

Moreover, clearly the Toda potential V(r) = e~ r + r — 1 satisfies the above assumptions. 

Theorem A. 2. Let x = (p,r) £ X and /i > 0. There exists a number v — v(fi, x) for which given any 
n, m € Z, the bound 



(A.5) 



Pn(t) 



r n {t) 



< q e -K\n-m\-v\t\) 



holds for all t £l and each z S {prmfm}, where 
(A.6) 

In fact, one may take 
(A.7) 



C = C(x) = max sup \V' (r n (t))\^ 2 , 1 

\(t,n)£lxZ . 



v = 2C[ e/ l+1 + - 



Proof. Fix x E X. Our previous theorem guarantees that for each x € M and n € Z, the function F„ : K — > M 2 
given by 



(A.. 



F„(i;x) 



r„(i) 



is well-defined and differentiable with respect to each z € {p m , r m }. When convenient, we will suppress the 
dependence of F n on x. Using the equations of motion, i.e. (IA.2I) . it is clear that 

(A.9) F n (t) = F„(0) + f ( , P "+\ (S) 7/ ( nis) . \ ds . 

Jo \ v { r n{s))-V'{r n -i{s))J 

Differentiating with respect to z we obtain 

(A.io) = a; F " (0) + 51 y o D n>e { s )—F n+e ( s )ds, 



with 
(A.11) 
Hence 
(A.12) 



D n ,e(s) 



|e|<l' 

5!(e)-5 (e) 
V'(r B )« (e) - V"(r n _i)*_i(e) 



oz 



< 



d 

7T^(0) 
oz 



I'l 



where 
(A.13) 



5i(e)+J (e) 
C 2 (5o(e) + 5_i(e)) 

Let us now consider the case that z = r m , i.e., 



(A.14) 



<9r 



' ? f„(o) (Tia 
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In this case, iteration yields 
d 



(A.15) 



dr, 



-F n (t) 



< 



< 



k=0 ' |ei|<l 



^ ] ^m+ejH 1 ■ i, 

|e*|<l 



Xn)D ei ---D ek 



E 

k—\n— m\ 



\2Ct\ k Dk (I 
kl 



where we have set 
(A.16) 



D 



ki<i 



c- 1 





Again, convergence is guaranteed since -§^F n {t) is continuous and thus bounded on compact time intervals. 
Using D 2 = H, one obtains 

\2Ct\ k 



(A.17) 



d 



dr, n 



F n (t) 



< C 



OC 



k—\n—7n\ 



The rest follows as in Theorem 12. II 

Note that in Flaschka variables 
(A.18) a n (t) = a n (t) (b n+1 (t) 



□ 



b n (t)) and b n (t) = -- (v'(- HK)) - v'(- HK-i))) 



the equations of motion read 
1 
2 

and (p, r) will be bounded if and only if (a, a" 1 , b) are bounded. The fixed point in these new coordinates is 
(do, bo) = (5,0) and (p, r) € X if and only if (a, &) — (ao, &o) = (a — | ; &) G 
Finally, by 

(A.19) a„(i) = a„(0)exp 

the sign of a n is preserved under the flow. 



(b n+1 (s) - b n (s))ds 
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